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Abstract. In this paper we apply to photoproduction total cross-section a model we 
have proposed for purely hadronic processes and which is based on QCD mini-jets and 
soft gluon re-summation. We compare the predictions of our model with the HERA 
data as well as with other models. For cosmic rays, our model predicts substantially 
higher cross-sections at TeV energies than models based on factorization but lower 
than models based on mini-jets alone, without soft gluons. We discuss the origin of 
this difference. 
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1 Introduction 

Cosmic ray experiments and planning for a fu- 
ture linear collider require knowledge of photon 
total cross-sections for values of very high cm. 
energy of the photon-proton or photon-photon 
system, in regions where no data are available. 
In this paper we study the predictions of a mini- 
jet model, extending to the photon a model de- 
veloped for proton-proton scattering [I] . We ap- 
ply this model to photo-production, leaving to 
a later paper the application to photon-photon 
scattering. 

Understanding the energy dependence of to- 
tal hadro-nic cross-sections continues to be an 
important issue in the study of strong interac- 
tions per se [2]. Over the years, various descrip- 
tions of this energy dependence have been given. 



Some approaches have focused on how far one 
can reach following the basic principles of analyt- 
icity, unitarity, factorisation etc., without any re- 
course to the details of the particular hadron in- 
volved, whereas at the other end of the spectrum, 
there are models which include the fundamentals 
of QCD as far as possible and then try to com- 
pute the cross-section in terms of measured prop- 
erties of the particular hadron. Of course, all de- 
scriptions have to be consistent with the require- 
ments of analyticity and unitarity. Most descrip- 
tions involve a few "soft" (non-perturbative) pa- 
rameters, which can not be determined through 
perturbative QCD. Again, basic symmetry, uni- 
tarity and factorisation arguments may at times 
lead to certain relationships among these soft pa- 
rameters for various hadrons. Often they may be 
determined only through fits to the experimen- 
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tal data and then one may only test approximate 
relations among these indicated by general argu- 
ments. In short, understanding the behaviour of 
the total hadronic cross-section and other soft 
quantities such as multiplicities etc., from first 
principles, is an extremely challenging problem 
and as stated before, one has different answers 
with varying degrees of relationship to QCD. 

Hadronic cross-sections for processes induced 
by the photon and the hadronic structure of the 
photon itself, have played a very interesting and 
important role, in furthering the attempts to un- 
derstand the theoretical issues involved in the 
subject. Photon-hadron interactions offer the the- 
orists one more laboratory to test their various 
ideas about computing "soft" quantities such as 
purely hadronic total cross-sections from basic 
principles. Historically, it is the interaction of 
the highly virtual photon with the hadron that 
offered the first glimpse of (almost free) quarks 
and later provided basic evidence for perturba- 
tive QCD being the correct dynamics to explain 
strong interactions in a certain kinematic do- 
main. However, in the present context, it is the 
photon structure function language i3j used to 
describe interactions of the real or quasi-real 
photon (invariant mass square ~ 0), with other 
hadrons or photon, that is of interest. In fact, 
the structure function of a quasi real photon at 
large values of x~f and that of a highly virtual 
photon (with large values of where —Q^ in- 
dicates the invariant mass square of the virtual 
photon) for all values of X'y , 00/11 be computed us- 
ing perturbative QED and QCD alone, for large 
values of momentum transfer square, of the 
probe. However, equally important is the (non 
perturbative) part of the real (or quasi real) pho- 
ton structure function at small which is not 
amenable to perturbative QCD (PQCD) compu- 
tations. 

In this paper, we apply our eikonal mini-jet 
model augmented by soft gluon resummation, 
which has been successful in providing an ac- 
ceptable description of the pp/pp data, to the 
description of total cross-sections of photon in- 
duced processes. In our model for the (purely 
hadronic) proton total cross-sections, we were 
able to compute the relevant components in terms 
of basic QCD inputs such as the experimentally 



measured parton densities and QCD subprocess 
cross-sections along with a few non-perturbative 
parameters. Given the prior success, it becomes 
of interest to see how the predictions of our model, 
applied to the total hadronic cross-sections of 
photon induced processes and using the experi- 
mentally determined knowledge on the structure 
of the "real" photon, compare with the data. We 
shall be mainly concerned with the issue of its 
energy dependence. 

To recapitulate: in this paper we explore the 
effects of the hadronic structure of the photon 
through studies of total cross-sections involving 
photons. While at low energy, these cross-sections 
can be obtained through factorization and vector 
meson dominance, the high energy range poses 
a different challenge. We have argued in a num- 
ber of papers [3,4,5,6] that the energy depen- 
dence of the photon induced processes do not 
follow from a straightforward application of fac- 
torization properties of the total cross-sections. 
We shall discuss various factorization results [H 
8,9,10,11] and compare some of them with the 
HERA data [T2l[T3] as well as with predictions 
of our QCD eikonal model with resummation, 
hereafter referred to as the BN model [14 . The 
reasons for this nomenclature will be clear as we 
describe the model. Some of its details are sum- 
marized in two Appendices, so as not to over- 
burden the reader with material published else- 
where. 



2 Total cross-sections: from pp to 77 

Experimentally, all total cross-sections rise asymp- 
totically with energy, but it is not yet clear whether 
the rate of increase is the same for different pro- 
cesses and whether their asymptotic behaviour 
is already controlled by the Froissart-Martin [TS] 
bound. For any given total hadronic cross-section, 
this bound says that asymptotically 

atot <C (log s)^. (1) 

From a theoretical point of view, this bound 
is only valid for the scattering of hadrons. At- 
tempts to extend it to virtual photon process 
[l6] have resulted in a less restrictive bound, 
but which cannot be extended to real photons. 
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leaving unanswered the question of whether the 
bound of Eq. ([1]) is valid also for real photons. 

Phenomenologically, the LEP data J_7j seem 
to indicate that the slope with which the total 77 
cross-section rises is not the same as in the pro- 
ton case 15^. This difference would spoil the sim- 
plicity of the so-called Regge-Pomeron model, in 
which the high energy rise is described through 
a single universal term [TT]. Of course, all total 
cross-sections do rise and to appreciate it at a 
glance, we show in Fig. [T]a compilation of data 
on pp/pp [18][19], 7p [l2lfT3l and 77 [17] scat- 
tering together with expectations from the BN 
model |14| for protons to be described in the next 
section. Since the data span an energy range of 
four orders of magnitude, with the cross-sections 
in the millibarn range for proton-proton, micro- 
barn range for photoproduction and nanobarns 
for photon-photon, to plot them all on the same 
scale, one needs a normalization factor. The data 
suggest to multiply the 7p cross-section by a fac- 
tor K, 330 and then 77 by (330)^, as shown in 

Fig. m 

It has been known for quite some time [20] 
that to get the photoproduction cross-section from 
the proton cross-sections in the region where they 
are approximately constant, namely after the ini- 
tial Regge-exchange type fall and before the be- 
ginning of the high energy rise, the multiplica- 
tive factor to apply for each photon leg in the 
cross-section can be obtained from Vector Me- 
son Dominance (VMD) (to go from a photon to 
a meson) and a quark counting factor, namely 
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With present p-meson data |21j and the relation 
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we would obtain i?-y « 1/360, consistent with 
the value indicated in the figure. 

Note that there is no a priori reason to expect 
the scaling factor to be energy independent. On 
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Fig. 1. Proton [19] and photon [T2][T3lfT7] normal- 
ized total cross-sections with the band expected from 
our BN model for pp/pp [I] and one typical curve 
from this model 14 . Cosmic ray data are from 
ref. [25]. 



the other hand, while at low energies the factor 
can be reliably evaluated through VMD, at 
high energies, it is likely to be different VF]\ due 
to the difference in the quark and gluon content 
of photons |3] versus that of the hadrons. The 
use of a simple multiplicative factor to compare 
the photon processes with each other and with 
the pure proton processes, is the simplest form 
of factorization. More complex forms of factor- 
ization exist in the literature, as in a recent for- 
mulation by Vereshkov and collaborators [53] or 
as in the model by Block et al. (also called the 
Aspen model) [24]. More comments shall follow 
in the last section. 



The above discussion points to the need for 
a description of high energy photon interactions 
where reliable predictions can be made based 
on the quark-parton structure of the photon. As 
stated earlier, we have developed such a model 
for purely hadronic processes [1, 14,26,27] which 
we shall extend and apply to photoproduction 
processes in the next section. 
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3 The Bloch-Nordsieck model (BN) 

This model is based on the eikonal representa- 
tion for the total cross-section [28„29j . and the 
eikonal incorporates QCD inputs such as parton- 
parton cross-sections, parton densities extracted 
from perturbative QCD fits to the data, actual 
kinematics, and soft gluon resummation. In de- 
tail, we use: 

1. QCD mini-jets [501151] to drive the rise of the 
total cross-section in the QCD asymptotic 
freedom regime; 

2. the eikonal representation for the total cross- 
section with the real part of the eikonal ap- 
proximated to zero and the imaginary part 
obtained through mini-jet QCD cross-sections; 

3. an impact parameter distribution obtained as 
the Fourier transform of the re-summed soft 
gluon transverse momentum distribution; 

4. resummation of soft gluon emission down to 
zero momentum to soften the rise due to the 
increasing number of gluon-gluon collisions 
between low-x, but still hard perturbative, 
gluons. 

While the eikonal representation with the mini- 
jet input has long been in use, our model differs 
from other existing eikonal models in that the 
impact parameter distribution is energy depen- 
dent and derived from soft gluon fcf-resummation, 
which gives the model its name. 

Before turning to the application of the model 3. 
to photon processes, we shall briefly discuss the 
various approximations involved. 

1. the energy dependence of QCD mini-jets re- 
flects the increase -at low x- of the gluon 
structure functions and this can provide a 
mechanism for the rise of total cross-sections. 
However, parton-parton cross-sections require 4. 
a minimum transverse momentum cut-off to 
avoid the divergence. Fixing such a cut- 
off around 1 GeV, leads to an extremely sharp 
rise. Our model starts with this input and 
then tempers the rise through the energy- 
dependent impact parameter distribution. 

2. the eikonal representation for total cross-sections 
is an approximation which allows one to en- 
force the requirement of s-channel unitarity. 



Here one obtains the total cross-sections through 
the eikonal formulae 

atot = '2. j (fb[l ^ e-'^"'^^''^'^ cosUex] (4) 
aei = j d26|i_e^x(M)|2 

a,^ei = j <fb[l - e-2a"x(b,.)] (6) 

The introduction of the jet cross-section as 
the term which drives the rise in the eikonal 
function can be done unambiguously through 
the inelastic cross-section, which only depends 
upon the imaginary part of the eikonal func- 
tion. Notice that the expression for ainei can 
also be obtained upon summing multiple col- 
lisions which are Poisson distributed with an 
average number n(b,s) = 2 SJm xib,s). Us- 
ing the experimental value of the p parameter 
(the ratio of the real to the imaginary part 
of the forward elastic amplitude) for pp and 
pp processes, we can estimate only about 4% 
correction from the real part. Hence, we as- 
sume SRex = 0, and thus obtain a very simple 
approximate expression 

atot = 2 J d^b[l - e-"('''^)/2] 

which can be used to test our ideas about the 
underlying strong interaction dynamics, 
the impact parameter distribution, which con- 
tains all the b-dependence, represents the mat- 
ter distribution in the colliding hadrons and 
our model obtains it as the shadow of the 
path followed by the partons as they scatter 
through the hadronic matter. The scattering 
process defines the partonic densities in the 
transverse plane. 

soft gluon emission is a QCD mechanism which 
introduces acoUinearity and thus can reduce 
the LO parton-parton cross-section. We use 
the resummation scheme obtained first in QED 
by Yennie, Frautschi and Suura [32] , extended 
to QCD transverse momentum (kt) distribu- 
tions in refs. [331134], applied to Drell-Yan 
pair and W- production in [35,36^, as well as 
to jet production in Earlier, we had dis- 
cussed such fct -resummation scheme for con- 
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stant but large values of the coupling con- 
stant in • As we discuss in Sect. 4, this 
resummation scheme includes an integration 
over the zero momentum modes. Depending 
on the (yet unknown) behaviour of as in this 
region, the integral over the low energy part 
of the single soft gluon distribution may be- 
come relevant. Similar integrals, involving as 
in the infrared region have been discussed 
in [39] as being related indeed to physical 
observables. For the transverse momentum 
case, this integral is usually absorbed in the 
intrinsic transverse momentum of the scat- 
tering hadrons and considered to be a con- 
stant. Our approach differs and has been dis- 
cussed in [1] and references therein. We shall 
return to this in Sect.H] 

The BN model was applied to proton-proton 
scattering, obtaining a total cross-section at LHC 
cr(Vi = 14 TeV) = 100 ± 12 mb, where the er- 
ror reflects various uncertainties such as in the 
choice of parton densities for the proton, mini- 
mum parton pt cut-off, called ptmin, and the in- 
frared behaviour of soft gluon coupling. Thus, 
the model has a number of parameters, some of 
which have a physical meaning associated with 
confinement. As such, we do not know how and 
if to change them as one goes from protons to 
photons. We shall try to vary them by no more 
than (5-10)% from their values for the proton. 
Whenever a stronger variation is required, we 
shall comment upon it. The model predictions 
are obviously dependent on the parton densi- 
ties in the photon: as in the case of the pro- 
ton, we have employed different available sets, 
obtained by fits to the data on the photon struc- 
ture function f^, and seen how best to describe 
the available data without appreciably changing 
the parameters. Application to photons however 
requires an additional insight. In the eikonal rep- 
resentation we need to adapt the hadronic lan- 
guage to that for the photon. One first needs the 
probability, Puad, that a photon behaves like a 
hadron and subsequently one may then use the 
eikonal representation, as in Refs. |241l401l¥T] : 

aZ = 2PnadJ d'hil-e--'"'^'-^/']. (8) 



Eq. [5] has no ab initio derivation. It is an ap- 
proximation, used in minijet inspired descrip- 
tions of photon cross-sections, useful to model 
the low energy term and provide a normalization 
for the cross-section. One could expect Phad to 
have an energy dependence or b-dependence, or 
both. However, to compare our model with other 
mini-jet models in the literature, and to make 
the simplest possible extension of our model from 
protons to photons, we consider it to be a con- 
stant as far as energy is concerned. As for the 
impact parameter space and related energy de- 
pendence, these effects are obtained in our model 
through QCD soft gluon resummation, just as in 
the case of protons. 

In Eq.[5]the real part of the eikonal has been 
approximated to zero, again following the proton 
model, while the imaginary part is obtained from 
the average number of inelastic collisions for a 
given impact parameter b, n'^P{b,s), at a given 
cm. energy ^/s. Following our BN model for pro- 
tons, we distinguish between collisions calculable 
as QCD mini-jets, and everything else, writing 
the average number of collisions as 

= <ftib, s) + A-^^ib, s)a]l,{s)lPuaa (9) 

with Ufiard including all outgoing parton pro- 
cesses with Pt > Ptmin- In Eq. [9]the impact pa- 
rameter dependence has been factored out, aver- 
aging over densities in a manner similar to what 
was done for the case of the proton in |27j. Be- 
cause the jet cross-sections are calculated using 
actual photon densities, which themselves give 
the probability of finding a given quark or gluon 
in a photon, Phad needs to be canceled out in 

'^hard- 

As for its value, Phad ~ Pvmd- Phad is 
not the same numerical factor used in Fig. [1] 
to normalize all the cross-sections at low energy, 
but it can be connected to it by making an ex- 
pansion of the eikonal in the low energy region, 
where a jet « 0, as shown at the end of this 
section. Also, while Phad can be factored out in 
some models, as we shall see later, this does not 
happen in the BN model. 

The mini-jet cross-section is obtained by in- 
tegrating the standard QCD inclusive jet cross- 
section, using a lower cutoff ptmin as described 
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Fig. 2. Left panel: Photon-proton jet cross-sections for different densities and a range of ptmin values. 
Right panel: average value of the maximum transverse momentum allowed for single initial state soft gluon 
emission, in 7p scattering. 



in Appendix A. The mini-jet cross-sections are 
to be calculated using parton densities (PDFs) 
for the proton and photon determined from per- 
turbative QCD analysis of the data on F2 as 
well as a variety of other data on hard processes 
for the proton. Common ones for the proton are 
GRV t42j, MRST CTEQ [44], whereas those 
for the photon are GRV[45], GRS [46], CJKL 
|47j . These densities are available both at lead- 
ing order (LO) or higher, but in our model we 
use only the LO ones, as part of the NLO ef- 
fects are described by soft gluon resummation 
and the use of NLO would result in some dou- 
ble counting. Of course, in using densities and 
parton-parton cross-sections only at LO but with 
resummation of soft gluons, our model lacks the 
non-infrared part of the NLO corrections. Since 
we consider the resummation effects in the in- 
frared region to be the most important for sat- 
uration and these are easily incorporated in our 
model, we have opted for LO densities, and thus 
also tree level parton-parton cross-sections and 
one loop as- We show in Figure O the energy de- 
pendence of the mini-jet cross-sections for 7p col- 
lisions, for two different sets of parton densities 
for the photon, GRS and CJKL. We have used 



different values of the cut-off, namely ptmin = 
1.2,1.3,1.4 GeV for GRS densities, higher val- 
ues for the case of CJKL densities, which give 
jet cross-sections rising faster with energy than 
those calculated using GRS [l^. As for the pro- 
ton densities, we have done all the model calcu- 
lations using GRV94. 

These cross-sections grow very rapidly with 
energy, reflecting the infinite range of QCD. Since 
the finiteness of strong interactions is reflected 
by the finite spatial extension of hadrons, one 
expects that the eikonal representation would 
check such growth through the impact param- 
eter distribution which appears in Eq. [51 A fre- 
quently used distribution is obtained as a convo- 
lution of the form factors of the colliding hadrons 
[49], namely 

A^Fib) = I -0^^%)^^{q)e"'''- (10) 

However, it was noted already in case of proton 
cross-section , that, without the inclusion of 
additional parameters, this choice is unable to 
reproduce both the early rise and the expected, 
Froissart -like, subsequent leveling off at high 
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energies. Apart from this purely phenomenologi- 
cal consideration, the form factor description be- 
comes undefined when dealing with photons. For 
photons, such models, labeled as Form Factor 
(FF) models, depend on how one defines the pho- 
ton form factor. In the literature, early attempts 
to apply the mini-jet eikonalized expression to 
the photon cross-sections [40] used a monopole 
expression for the photon (as in the pion case) 
and the usual dipole expression for the proton 
form factor with = 0.71 GeV^, obtaining 



1 



Ait kl 
kl-v 



[vbKi (vb) 



rilCoii^b) - ICoikob))] (11) 



with /cq = 0.44 GeV^ . The above expression 
can be adapted to photon data by varying the 
parameter fco, and then the pion form factor ex- 
pression for the photon can be understood to 
represent an intrinsic transverse momentum |50t 
[5T| . In the Aspen model [H] another possibil- 
ity has been explored. Namely, the overlap func- 
tion is parametrized as the Fourier transform of 
a dipole form factor 



w{b,^^) = ^{^,b)'K,{^,b) 

9D7r 



(12) 



with three different scaling parameters for the 
three terms in which the eikonal is split, quark- 
quark, quark-gluon or gluon-gluon scattering. In 
this model, one then uses a single functional ex- 
pression for the b-distributions in hadron-hadron, 
hadron-photon or photon-photon scattering, but 
the difference between these different processes 
is provided through the fi parameters which scale 
among the various processes according to the 
additive quark model. In another QCD inspired 
model [52]. a similar modelling has been made. 

More fundamental attempts to obtain the pho- 
ton impact factor in the context of perturbative 
QCD can be found in 53] and references therein. 
However, since these models are derived for vir- 
tual photon scattering, they cannot be used in 
the context of our BN mini-jet model, which 
is applied to real photon processes, with real 
photon parton densities. In this paper we fol- 
low the same strategy which we used in case of 



the proton cross-sections and, for hard collisions, 
use mini-jets and soft gluon resummation, with 
nhard given by: 



,7P 



Al%{b,s)a''^ 



jet 



Ph 



iad 



(13) 



with 



2^ d2p(K^)^_,K,.b 



e 



-h(b,qmaa:) 



= ABNib,qmax{s)). (14) 



The function A^j^ is normalized to 1 and is ob- 
tained from the Fourier transform of the soft 
gluon resummed transverse momentum distribu- 
tion, whose structure we discuss in the next sec- 
tion. 

To complete the calculation of nj^rd^ '^^^ 
to specify the value of Phad, which in eikonal 
models [24ll4T] indicates the probability that a 
photon behaves like a hadron and is defined by 
the low energy part of the cross-section. At low 
energy, namely for ^/s « 5 10 GeV, the mini- 
jet cross-section is very small and n{b, s)^^ ~ 



'soft 



(&, s). This part of the cross-section is out- 
side the range of perturbative QCD we have de- 
scribed so far. Using Eq. [51 we find that we can 
get a good description of the low energy data 
for the total cross-section with 



•'soft 



2 



^pp 



(6,s) 



(15) 



where JT-so/tl^i same function we have 

used for our description of proton-proton colli- 
sion in ref. [T] and Phad = 1/240, a result con- 
sistent with Eq. [5] 



4 The impact parameter distribution 
and the saturation parameters 

The distribution Abn is energy dependent through 
the quantity Qmaxis), which represents the aver- 
age maximum transverse momentum allowed to 
a single soft gluon emitted in the initial state in 
a given hadronic collision. This quantity is the 
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input to the kernel h{b, qmax), which describes 
the exponentiated, infrared safe, number of sin- 
gle soft gluons of all allowed momenta and is 
given by. 



Hb, qmaxis)) = 

X (lot 



16 dkt as{k^) 

2q 

max is) 



[l-Mhb)] (16) 



In the BN model, this function provides the cut- 
off in impact parameter space which softens the 
rapid rise of the mini-jet cross-section. To render 
this point clear, we shall summarize and outline 
our argument in what follows. 

The soft gluon resummation formula in the 
transverse momentum variable has been known 
for a long time and reads [33ll34l[38] : 



with 



(17) 



M6,w(s))= T'"""'^ ^d'fi(fc)[l-e^''-b] (18) 
Jo 

where Qmaxis) is the maximum transverse mo- 
mentum kinematically allowed for single emis- 
sion. The expression in Eq. [T7] can be obtained 
from the more general expression for soft resum- 
mation 



d*P{K) = d'^K 



d^ 



JK-x~h{x,E) 



(19) 



with h{x, E) similarly defined as in Eq. [TH] This 
expression was obtained long time ago in QED 
[32] J with order by order cancellation of the in- 
frared divergence in perturbation theory. Notice 
that the same expression is obtained very sim- 
ply in a semi-classical way using the methods of 
statistical mechanics and imposing energy mo- 
mentum conservation [54j . In such derivation, it 
is energy-momentum conservation that brings in 
the cancellation of the infrared divergence be- 
tween soft and real quanta emission. 

In QED, d^n{k) oc alog( ^^^'";" ) and re- 
summation in the transverse momentum vari- 



able is well approximated by a first order expan- 
sion in a. For large values of the coupling con- 
stant however, this approximation would be in- 
adequate and we noticed long time ago [38] that, 
for such cases, resummation of soft quanta emis- 
sion can provide a transverse momentum cutoff, 
based on the expansion of the Bessel function 
Jo{ktb) around fc(6 w 0. 

In QCD [33,34J, Eq.[TH| shows the presence of 
a non-trivial complication, namely the impossi- 
bility to use the asymptotic freedom expression 
for as down to fcj = 0. It is this non-trivial com- 
plication which we exploit to study scattering in 
the very large impact parameter region. 

To overcome the difficulty arising from the 
infrared region, the function h(b,E), which de- 
scribes the relative transverse momentum distri- 
bution induced by soft gluon emission from a 
pair of, initially coUinear, colliding partons at 
LO, is split into 



h{b, E) = co(m, b, E) + Ah[b, E), 



(20) 



where 



Ah{b,E) 



16 /-^ asjkt) ,, ji,,..dh 2E 

T / [l-Jo(6fcf)]— In— . 

3 7^ TT kt kt 

(21) 

Since, in Ah(b,E), the integration only extends 
down to the scale fi (not zero), the Jo{bkt) is as- 
sumed to oscillate to zero and hence is dropped. 
The last integral is now independent of b and can 
be performed, giving 



Ah{b,E) 



32 



33 



ln(. 



,2E 
A ' 



ln(ln(|))-ln(ln(^)) 



2Nf 
■ (22) 



where A being the scale in the asymptotic free- 
dom expression for as- In the range 1/ E < b < 
1/ A the effective h^f / (6, E) is obtained by set- 
ting fjL = 1/b [31]. This choice of the scale intro- 
duces a cut-off in impact parameter space which 
is stronger than any power, since the radiation 
function is now 



f(b,E) _ rMlAjl^](16/25);n(£;V^') 



-ln{E^/A^ 



(23) 
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which is Equation(3.6) of ref. [31]. The remain- 
ing b-dependent terms in E) are dropped, 
a reasonable approximation if one assumes that 
there is no physical singularity in the range of 
integration < kt < 1/b. However, when the 
integration in impact parameter space extends 
to very large b- values, as is the case for the cal- 
culation of total cross-sections, this expression 
fails to reproduce the entire range of the energy 
dependence of low energy transverse momentum 
effects. To explore the very large b-region, we 
suggest to use Eq. with its full integration 
range, proposing a phcnomcnological approach 
to the zero-momentum soft gluons 55 . Our ap- 
proach uses a singular, but integrable expression 
for as- this allows us to extend the integral to 
the minimum allowed value zero and to obtain a 
b-dependence which is stronger than the one of 

Thus the infrared region can provide an im- 
pact parameter cutoff at large b- values, provided 
the integral is finite. This requires any proposed 
expression for as in the infrared region to be in- 
tegrable. As discussed in detail in ref. [IB], the 
actual functional dependence of the cut-off de- 
pends on the model for as in the infrared. Here 
we mention that in the model we propose, the 
resulting cut off in b-space is at least an expo- 
nential function. 

In the next subsections, we shall discuss var- 
ious proposals to model the infrared region and 
how the full integral of Eq. (fT6| controls the 
saturation of the cross-section through its limits 
of integration. We shall see that our model for 
soft gluon emission is regulated by a constant 
infrared parameter p and the energy dependent 
momentum function Qmax as follows: 



1. the energy dependent momentum saturation 
parameter qmax{s) depends on the energy be- 
haviour of the density functions of colliding 
partons and on ptmim the mini-jet cut-off, 

2. the infrared parameter p, to be specified short- 
ly, defines the infrared behaviour of asik^). 
The closer its value is to 1, the more the 
mini-jet cross-sections will be quenched at 
any given energy. 



4.1 The momentum saturation parameter 

Qmax (-5) 

For any given parton parton collision, qmax{s) 
can be defined by kinematics. We introduced this 
quantity for the first time in |27j to represent 
the maximum transverse momentum carried by 
a single gluon, averaged over the basic scattering 
cross-section with a procedure described in Ap- 
pendix B for the convenience of the reader. To 
highlight the physical meaning of qmax{s), let us 
define the saturation parameter 

for each parton pair of cm. sub energy s which 
scatters into a final parton pair of cm. energy 
■y/ Sjets- Let us now use the kinematics of the 
process 

parton{xi)+parton(x2) — > gluon{kt)+jeti+jet2 

(25) 

to write the maximum transverse momentum of 
the emitted gluon, in the case of limited energy 
loss as [55] 

ktmax ^ ^{l- ^ (26) 

2 s 2 

This quantity plays an important role in our 
model. As the available cm. energy increases, 
it increases, depending upon the probability of 
producing a parton pair scattering into a given fi- 
nal state. It thus depends upon the densities and 
the parton-parton cross-section. As it increases, 
more and more acollinearity is introduced in the 
scattering and the stronger is then the reduction 
in the growth of the mini-jet cross-section. 

Notice that now there appear two different 
scales and both low-x perturbative gluons as well 
as soft gluons. We stress the distinction between 
them: low-x gluons participate in the hard parton- 
parton scattering described by the mini-jet cross- 
section discussed in the previous section, for which 

Ptout = > Ptmin w 1 4- 2 GeV (27) 

These low-x perturbative gluons interact with 
a strength proportional to ctsiPtout)^ while soft 
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gluons are those emitted, from the initial state, 
in any given parton-parton process with trans- 
verse momentum 



kt < kt 



10 ^ 20% ptout 



(28) 



This scale, ktmax defines the single soft gluons, 
whose number can be indefinite. These soft glu- 
ons need to be re-summed through the proce- 
dure which results in the exponentiated factor 
of Eq. [H 

In a model such as ours, which is not a Monte 
Carlo simulation of the processes involved, we 
have opted for averaging these effects, embody- 
ing them in a factorized expression such as that 
given by Eq. [141 with ktmax averaged out to ob- 
tain qmax, as shown in Appendix B. The ex- 
pression for qmax(s) depends both on the par- 
ton densities and the value of ptmin- The result- 
ing quantity is energy dependent since the den- 
sities are energy dependent through the applied 
DGLAP evolution. The averaging process done 
in this model includes only quark densities as 
the source of the leading acoUinearity effect. We 
consider the leading effect to arise because of 
soft gluon emission from the external legs of the 
scattering process, valence quarks for the proton 
beam and all flavours of quarks for the photons. 
An improvement of the model could include soft 
gluon emission also from the low-x perturbative 
gluons, as we shall discuss in a forthcoming pa- 
per. In the right-hand panel of Fig.[2]we show the 
dependence of Qmaxis) upon the cm. energy of 
the colliding particles, for the same densities and 
Ptmin values used in the mini-jet cross-sections 
shown in the left panel. 

As Qmax incrcascs with energy, the growth of 
the total cross-section due to mini jets is tem- 
pered by soft gluon emission, through the ex- 
ponential damping factor e-Hb,qma^) _ However, 
there is an equilibrium between the increase of 
Qmax and the rate of increase of the mini-jet 
cross-section since one reflects the quarks and 
the other the gluon densities. The distribution 
of these partons at high energy follows the par- 
ton sum rules and one is not independent of the 
other. From the right hand panel of Fig. [2] we 
see that Qmax, for both GRS and CJKL densi- 
ties, will reach some sort of saturation at high 



energies, which reflects in the total cross-sections 
reaching a stable slope. 

The momentum saturation parameter qmax 
is not the only quantity which gives rise to satu- 
ration, the infrared limit of also plays a major 
role. We shall discuss this in the next subsection. 



4.2 A phenomenological approach to the 
infrared limit of as 

To complete the calculation of the impact pa- 
rameter distribution for hard processes in jp col- 
lisions, we need to discuss the lower limit of inte- 
gration in Eq[T6l Usually, the soft gluon resum- 
mation formula extends the soft gluon momenta 
to an infrared cut-off taken to correspond to the 
intrinsic transverse momentum scale of the scat- 
tering hadrons I33II34! . Instead, in our model, we 
extend the integration down to the zero momen- 
tum modes. 

The reason to do so lies in the nature of the 
cancellation of the infrared divergence. To obtain 
this cancellation, it is mandatory that virtual 
and real soft gluon emission join in the zero mo- 
mentum limit. Of course we do not know how to 
deal with as in this limit and for this reason the 
integral over this region is usually left out and 
substituted with a constant intrinsic transverse 
momentum. However, this part of the integral 
is relevant for many minimum bias processes, as 
pointed out in ref. [39 . With our model, we aim 
to relate the behaviour of the impact factor at 
very large b-values with the infrared region of 
soft gluon emission. 

To do so, we need therefore to make an ansatz 
as to the behaviour of the strong coupling con- 
stant in the infrared region, where the usual asymp- 
totic freedom expression for as(Q^) cannot be 
used. One possibility is to use an expression which 
would go to a constant as — > as in 



as{Q^) 



127r 



1 



33-2A^/log[a+^] 



(29) 



with a w 2 [57(135] and A = Aqc d ■ This expres- 
sion is often referred to as the frozen as case. 
Another possibility is to employ the Richardson 
potential for quarkonium, which uses a singular 
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as, namely a = 1, so that 



(30) 



The Richardson potential has been shown to give 
good results to describe charmonium states [SS] , 
but it cannot be used here because the integral 
over the soft gluon modes would diverge. The 
reason it works in quarkonium applications is 
that in that case one never actually reaches val- 
ues corresponding to = 0, since the potential 
binds the two quarks in a region of space at fixed 
finite distance of 0{rBohr)- 

In order to be able to use the Richardson-like 
af', we soften the singularity with the proposal 
that in the infrared limit, one can phenomcno- 
logically use the expression 



Cis{kt) — constant ^ ( ^ ) 



fct ^ (31) 



where yl is a cut-off of order Aqcd, and p is a 
parameter which embodies the infrared behav- 
ior, with p < 1 so that the soft gluon integrals 
converge. For the time being, we consider the 
above expression as a phenomenological ansatz . 
The constant in front of Eq.[3T] should be chosen 
to provide a smooth extrapolation to the per- 
turbative expression for a^. Our choice for the 
interpolating function is 



127r 



a. 



33 - 2Nf ln[l 



(32) 



This expression was also introduced to describe 
the intrinsic transverse momentum of Drell-Yan 
pairs, with the choice A = 100 MeV [55] and 
p = 5/6. This choice for the infrared behaviour 
(zero momentum gluons) was motivated [14j by 
an argument due to Polyakov [59] . It is clear that 
the closer p is to 1, the bigger the soft gluon in- 
tegral h{b,qmax{s)) is and the stronger the sat- 
uration effects will be. 

We shall show the results for the total 7p 
cross section for this and other models in the 
next section. 



5 Total 7p cross-section at accelerator 
energies 

In this section we examine the data on cr*p*, 
starting from the low energy photoproduction 
data [ISj up to the high energy HERA data and 
compare them with model predictions. In fact, 
HERA studied ep scattering and not only yielded 
information on energy dependence of the photo- 
production cross-sections but also on 7*p cross- 
sections for different values of virtuality (invari- 
ant mass —Q^) of the photon. For small values of 
cc-y and Q^, the measurement of F2{W'^ , Q^)/Q'^, 
where W is the invariant mass of the 7*p sys- 
tem, gives a measurement of a"^ ^ through the 
relation, 

The extrapolation of the measured quantity on 
the right hand side can then give the photo- 
production cross-section in the limit — 
[50] . The extrapolation of data collected with the 
ZEUS Beam Pipe Calorimeter (BPC) [HI], based 
on a Generalized Vector Meson Dominance model, 
produces |62p63j a set of measurements in a con- 
tinuous energy range W^p = 104 -^251 GeV 
consistent within errors with the photoproduc- 
tion data. The systematic errors are due to the 
GVMD extrapolation, with the longitudinal cross- 
section (Ti = 0, of the ZEUS BPC95 data. In the 
comparison, we have also included some cosmic 
ray data [64j which are in an energy range lower 
than HERA data, but higher than the older pho- 
toproduction data. 

Let us start with the BN model for photons 
as described in the previous section. We have 
used GRV densities for the protons 42J and have 
varied the photon densities, using both GRS and 
CJKL. We show the result of the model and the 
dependence upon the model parameters in Figs. 
I3I4I In FiglSjwe have varied ptmin and the densi- 
ties to describe the high energy data from HERA 
in addition to the most acceptable description of 
the beginning of the rise, while keeping the pa- 
rameter p in a range close to the pp/pp case. In 
Fig. |4| we have allowed for a larger variation in 
the value of the infrared parameter p, fixing the 
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PDF set and a range of appropriate values for 

Ptmin • 

In order to obtain a good model description, 
we shall focus not only on the HERA data, but 
also on the beginning of the rise, as this signals 
the onset of the contribution of QCD processes 
and is strongly dependent upon ptmin- We can 
see from Fig.[2]that, for the range otptmin values 
of interest, the mini-jet cross-sections calculated 
with CJKL densities rise faster than those cal- 
culated with GRS. It follows that, to describe 
the same HERA data, one will need to use dif- 
ferent values of ptmin depending upon the PDF 
set used. Thus CJKL densities call for a larger 
Ptmin than GRS densities. In Fig. [3] the infrared 
parameter p has been kept close to the value de- 
termined from the pp/pp cross-section, namely 
p « 0.7 0.8. We see that the range of ac- 
ceptable Ptmin values for GRS densities is not 
far from those used in the pp/pp case, where 
Ptmin ~ 1-1 ^ 1-25 GeV, but it is higher for 
CJKL. 

To summarize the results of these figures, 
the latest HERA data are well described for a 
range of parameters p = 0.75 0.8 and ptmin ~ 
1.2 1.3 GeV, to be compared with the pp and 
pp case where the range was very similar, with 
our central value p = 0.75 and ptmin = 1-15 GeV 
for GRV densities. A good description is also ob- 
tained with CJKL densities, but then one needs 
a different {p-,ptmin} set, as one reads from the 
second panel in Fig. [S] 



5.1 Parameter dependence 

In the previous figures, we have applied the BN 
model for protons to jp scattering. No attempt 
has been made to choose one particular curve 
through an evaluation of the of our curves 
with respect to the data. This is so because the 
aim of this paper is to extend the BN model to 
photon processes with as little changes as possi- 
ble from the proton case. Another reason is that 
the data have rather large errors, and are ex- 
tracted from very different experimental situa- 
tions, which range from cosmic rays to HERA 
data on photoproduction to extrapolations of data 
cohected with the ZEUS BPC. We did not wish 
to select one particular data set and adopted the 




Photoproduction data before HERA 
ZEUS 96 
HI 94 

Vereshkov 03 
ZEUS BPC 95 



GRS(7) GRV94(p) with BN resummation 
Soft from pp soft eil<onal times 2/3 

— p=0.75andp,„i„=1.2,1.3,1.4GeV 

— p=0.8 and p„i„=1 .2,1.3 GeV 



^° Vs(GeV) 




Photoproduction data before HERA 
■ ZEUS 96 
• HI 94 

i Vereshl<ov 03 
□ ZEUS BPC 95 



CJKL(y) GRV94(p) with BN resummation 
Soft from pp soft eil<onai times 2/3 

— p=0.75andp,^i„=1.6,1.7,1.8GeV 
p=0.8 and p,„i„=1 .6,1.7,1.8 GeV 

^° Vs(GeV) 



Fig. 3. Total 7p cross-section with a range of param- 
eter values close to the proton case, GRV densities 
for the proton and GRS or CJKL densities for the 
photon. Data from HERA are from Zeus [T3 ] . HI p ] 
and a set of data from the ZEUS BPC extrapolated 

from qv [oniiniES]- 
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strategy to see whether the range of parameter 
values used for the proton results give curves 
consistent with the photo-production data. 

The high energy behaviour of the model de- 
pends on the following entries: 

1. PDF set, for which we have used GRV for the 
proton and GRS or CJKL for the photon 

2. minimum hard parton cutoff ptmin 

3. infrared parameter p 

with the choice of ptmin related to the chosen 
PDF set. Of the two parameters, ptmin is ba- 
sically fixed so as to reproduce the early rise, 
where soft gluon resummation is not yet impor- 
tant, while the parameter p controls the quench- 
ing of the rise at high energy and also the ab- 
solute value of nfiard{b, s). We also use ptmin ~ 
1-^2 GeV for the perturbative mini-jet calcula- 
tion to make sense. The dependence on densities 
and Ptmin was shown in Fig. [31 while the parame- 
ter p was kept within the range of values used for 
the proton case. We see that the GRV/GRS den- 
sity set is the one for which the extrapolation of 
the model from protons to photons works better. 
On the other hand, the model is also able to de- 
scribe the data using CJKL densities, but in such 
case, the parameters are different from the pro- 
ton scattering case. In Fig. |4] we give examples 
of how it is possible to quench the high energy 
rise from a given choice of PDFs (CJKL, in this 
example) through the parameter p. Notice that 
to catch the early rise, around ^/s — 20 GeV, 
for CJKL densities one needs ptTnin ~ 1-5 GeV, 
but then this requires a larger p-value in order 
to quench the rise and not overshoot the HERA 
data points. 

All in all, we can say that the model ade- 
quately describes the photon-proton cross-section 
data and we can try to extend it to higher ener- 
gies so as to make predictions for cosmic ray en- 
ergies to be reached by the AUGER experiment 
[65,66J. We turn to this problem in the next sec- 
tion. But before this, we address the question of 
factorization: is a photon like a proton just mul- 
tiplied by a constant factor? From what we have 
seen so far, one could describe 7p total cross- 
section up to HERA energies either through a 
microscopic model such as our BN model, with 
quarks and gluons, or through other approaches 
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l°Vs(GeV) 10 

Fig. 4. Total 7P cross-section with GRV for proton 
and CJKL densities for the photon, for a spread of 
p values. 



based on various forms of factorization. In par- 
ticular, the Aspen model also gives a good de- 
scription as do other approaches, based on mul- 
tiplying the result of fitting pp/pp data with a 
constant factor. We shall discuss this point in 
the coming subsection. 



5.2 Factorization: a hadron-like photon 

In the previous section, we have applied our model 
to the total jp cross-section, using available pho- 
ton densities, going through the various steps 
defining our model, namely calculation of mini- 
jet cross-sections, evaluation of the energy de- 
pendence saturation parameters, determination 
of the energy dependent impact parameter func- 
tion from soft gluon resummation AsNib, qmax{s)) 
and finally eikonalization. In this approach, at 
high energy, the photon is an independent en- 
tity from a hadron, with the rising behaviour of 
the cross-section and the b-distribution of the jp 
collision determined independently from other 
hadron — hadron collisions such as pp. This is 
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different from other models, for instance the As- 
pen model |24) , where the photon properties are 
obtained through scaling factors inspired by the 
additive quark model. As a consequence, in the 
Aspen model for photons, one can prove a fac- 
torization property 8, which would then allow 
to extract the 77 cross-section simply as [7] 



(33) 



with cr„ to indicate the nucleon cross-sections. 
We shall discuss the 77 cross-sections within our 
BN model for photons in a separate paper, how- 
ever we notice that such factorization is not to 
be expected in the model we present here. 

Other types of factorization models are based 
on the Regge-Pomeron exchange, keeping a con- 
stant universal behaviour of the rising part of 
the cross-section with coefficients based on the 
factorization of the residues at the poles in the 
elastic amplitude, so that 



-yn 



7n^ 



_77 _ 



Y 



(34) 
(35) 

(36) 



with e K, 0.08 0.09. This type of factorization 
is of course different from the one in Eq[331 but 
it still implies the idea that there is a universal 
behaviour of the energy dependence, not only 
at low energy, where one can confidently assume 
that the hadronic interactions of the photons are 
those of a vector meson, but also at high energy. 

Such a description of the photon, i.e, that the 
photon is always hadron-like, could be reflected 
in our model by simply scaling the BN cross- 
section for protons, as 



"tot — ^l^tot 



d'^b[l 



-n''P(6,s)/2i 



(37) 

Present accelerator data for jp are consistent 
with factorization models, including an applica- 
tion as given in Eq. I37[ but as we shall see in the 
next section, at higher energies, expectations will 
differ. 



6 Extrapolation to very high energies 



In this section we extend our calculation be- 
yond present accelerator energies and compare 
our predictions with other approaches. We start 
with the simplest factorization model of Eq. [37] 
and multiply the band of results obtained in ref. 
[Tl for proton-proton total cross-section with a 
constant factor. This is similar to what we did 
in Fig. [U except that we use the full band from 
Fig. 2 of ref. ^1] . Let us indicate these predictions 
as BNp = B N protons /•ii^ (F for factorization). 
We then compare this band with the results ob- 
tained using the BN model with photon densi- 
ties, GRS and CJKL, namely the curves shown 
in Fig. H extended to ^/s:;^ = 20 TeV. This 
comparison is shown in Fig. [51 We see that, at 
energies around and through the TeV region, the 
band obtained from af^j falls short of what the 
BN model for photons {BN^) predicts. Other 
models, which enjoy factorization like the Aspen 
model, also remain lower than our curves. While 
at moderate, HERA like energies, all the three 
models, Aspen, BNj or BNp give acceptable fits 
to the data, there is a difference of almost 50% 
among their high energy extrapolations. Thus, 
the first interesting conclusion from this explo- 
ration of the very high energy region is that there 
is a distinct difference between predictions from 
our BN model and those from the QCD inspired 
model of Block et al. (Aspen) [M] , as well as from 
a straightforward multiplication of our band of 
predictions for the proton times a normalization 
factor. 

The next interesting result from this extrap- 
olation appears when one compares our model 
predictions with the fit to HERA data performed 
by Block and Halzen and based on a low energy 
parametrization of 7p resonances joined with Fi- 
nite Energy Sum Rules (FESR) and asymptotic 
In'^ s behaviour 67]. Fig. [Sj shows a band cor- 
responding to the predictions of our model for 
photons (upper band) compared to BNp (lower 
band), the Block and Halzen fit [67j, the As- 
pen model of [241, and an eikonal mini-jet curve 
which uses the proton and pion form factors for 
the impact parameter distribution (FF model). 
The central (full ) curve in the upper band corre- 
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Fig. 5. Total 7p cross-section with GRS and CJKL 
densities, compared with ref. [24] predictions and 
with a brute force factorization of our proton-proton 
results from pj. 
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Fig. 6. Total 7p cross-section with expectation from 
BN model using GRS and CJKL densities (upper 
band) for the photon, compared with the model in 
24 , the fit described in [67], a factorization model 
(lower band) and the eikonal mini-jet model without 
soft gluons (dot-dashes). 



Fig. [6] deserves some comment. For the curves 
shown in this figure, the parameters have been 
chosen so as to reproduce the highest available 
accelerator data (through ptmin and p values for 
the BN model, and through ptmin for the FF 
model) and the low energy data, the latter through 
Phad and ctq. As the cm. energy increases, the 
model results show noticeable differences between 
the hadron-like models, Aspen and BNp, and 
the photon-density model BN-y, and much more 
between all of them and the eikonal mini-jet (EMM) 
Form Factor model. Neglecting the FF model, 
which we think is incomplete, we nonetheless 
have a remarkable difference in the very high en- 
ergy range, 10 TeV and beyond. Because these 
predictions may impact strongly on the photon 
content of high energy cosmic rays [65il66j . this 
difference does matter. 
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We notice that the curve, labelled Block-Halzer 
(BH), from [67] lies within the band of the BN-y 
model. The BH curve is based on a best fit to 
low energy 7p data, joined smoothly with a fit 
of high energy accelerator [T2jfT3] and cosmic ray 
data [64 of the form 

cTyp = co+ci log(t//m)+C2 log^{iy/m)+(3-p>/y^v/m 

(38) 

where v s the laboratory photon energy. There is 
a noticeable difference between the slope in the 
rising part of the cross-section between the As- 
pen model and the BH fit, as there is between the 
modelling content between all these descriptions. 
In our model the rise is based on the gluon den- 
sities entering the calculation of the QCD mini- 
jets cross-sections and on the soft gluon resum- 
mation ansatz for the impact parameter distri- 
bution. The calculation of these inputs relies on 
realistic PDF distributions and actual, LO, par- 
ton parton cross-section. Then, the very high en- 
ergy (in the TeV region) agreement between the 
BH best fit based on an analytic expression and 
our results is an independent check of the correct 
physics content of the BN model. This fit con- 
firms the inherent interest of our approach based 
on QCD mini-jets and soft gluon resummation. 

For possible use, we report in table [U the 
numerical values obtained in our model for the 
cross-sections shown in Fig. [51 

From a numerical point of view, the curves 
for the BN model for protons and photons differ 
because the b-distributions for protons or pho- 
tons obtained in our model from ABNi,qmax{s), b) 
differ, and they differ because the maximum mo- 
mentum allowed to individual soft gluons is dif- 
ferent. This quantity for the hard part is ob- 
tained through the kinematic constraint aver- 
aged over the quark densities and the latter are 
of course different for protons and photons. This 
is apparent from a comparison of qmax{s) for 
pp [UlSSj and 7p : for comparable cm. energies 
qmax{s) for pp rises to higher values than the one 
for 7p, resulting in more saturation for pp. These 
differences are due to the quark densities enter- 
ing the averaging process defining qmax{s): den- 
sities are a phenomenologically extracted quan- 
tity and as such it is to be expected that they 
reflect the different structure of the interacting 
particles, namely the difference between valence 



quarks bound in a proton and quark pairs in 
which the photon will split and their respective 
evolution. 



7 Conclusions 

We have applied to 7p scattering an eikonal mini- 
jet model with soft gluon resummation devel- 
oped for the proton total cross-section. The model 
relies on the parton structure of protons and 
photons and indicates a different high energy be- 
haviour for 7p relative to pp and pp. We sug- 
gest that this different behaviour may be due 
to the different parton structure and high en- 
ergy evolution properties of quarks in the proton 
and quarks in the photon. Furthermore, this re- 
sult strengthens our confidence in the BN model 
as a good approximation to a QCD description 
of hadronic interactions in minimum bias pro- 
cesses. 
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Appendix A: The mini-jet cross-section 

The QCD jet cross-section for the process 

hadroHA + hadrons — » X + jet (Al) 

is obtained by embedding the parton-parton sub- 
process cross-section with the given parton den- 
sities and integrating over all values of incoming 
parton momenta and outgoing parton transverse 
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Table 1. Values (in mb) for total cross-section for "fp scattering evaluated in the cm. energy of colliding 
particles, corresponding to the bands shown in Fig. 



\/5 

GeV 


EMM with Form 

Factors,GRS 
Ptmin = 1.5 GeV 


BN-y model 
(upper curve) 
top band 


BN-y model 
(lower curve) 
top band 


B Nproton/SSO 

(upper curve) 
lower band 


B Nproton/SSO 

(lower curve) 
lower band 


5 


0.116 


0.116 


0.116 


0.118 


0.119 


10 
11.46 


0.114 


0.115 


0.114 


0.115 


0.116 


48.93 
50 


0.122 


0.130 


0.121 


0.131 


0.129 


100 
112.14 


0.139 


0.155 


0.140 


0.15 


0.143 


478.74 
500 


0.238 


0.228 


0.203 


0.199 


0.182 


1000 
1097.3 


0.352 


0.279 


0.250 


0.221 


0.199 


4684.6 
5000 


0.635 


0.384 


0.338 


0.280 


0.240 


9000 








0.310 


0.255 


10736.8 


0.829 


0.449 


0.390 






14000 








0.335 


0.266 


20000 


0.985 


0.499 


0.429 







Appendix B: The calculation of qmax{s) 

Simple kinematics can give the maximum trans- 
verse momentum allowed to single giuon emis- 
sion in a process like 

partoni{xi) + partcm,2{x2) — > gluon{k) -\- X{Q) 

(Bl) 

namely 

M(xi,x2,Q2)-^(l-?) (B2) 
2 s 

with s = sx\Xi- If X represents two jets from 
the outgoing parton-antiparton pair, one can use 
Q2 k, 4pj . The calculation is simplified by intro- 
ducing an average over the parton parton cross- 
section and integrate over all x values [SS] ob- 
taining 



momentum pt , according to the expression 



a-^^{s, Ptmin) = dpt dxi dX2 

dd^Hs) 

X hAixuP^t)hBix2,P^t)^^ (A2) 

i,j,k,l ^* 



where A and B are the colliding hadrons or pho- 
tons, in this case A— proton, _B — 7. By construc- 
tion, this cross-section depends on the particu- 
lar paramctrization of the DGLAP [53] evoluted 
parton densities, some of which do extend to very 
low x-values but not too high p^ values. This 
cross-section strongly depends on the lowest pt 
value on which one integrates. The term mini-jet 
was introduced long ago [301I3T] to indicate all 
those low Pt processes which one can still expect 
to be QCD calculable but which are actually not 
observed as hard jets, pt being the scale at which 
to evaluate as in the mini-jet cross-section cal- 
culation, one can have ptmin ~ 1 2 GeV . 



/ [dxidx2) J^jdzy/xiX2{l - z)D{xi,X2) 



^ J {dxidx2) dzD{xi,X2) 

(B3) 
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where 2„im — ^Ptmin/^^ ^ denotes the usual 
quark density expression 



11. 



A. Donnachie and P. V. Landshoff, Phys. 
Lett. B 296 (1992) 227 arXiv:hep-ph/9209205; ; 
A. Donnachie and P. V. Landshoff, Phys. Lett. 



Dix,,X2) - y[Mxi)/x,][f,{x2)/x2] (B4) B 595(2004) 393 [arXiv:hep-ph/0402081 



and we have also assumed that the parton-parton 
cross-section, appearing at both numerator and 
denominator, can be evaluated at its maximum 
value, pt — Ptmin, thus dropping out of the cal- 
culation. 
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